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SUI....3VIARY 


This report eontaras a description for a sinrulation program for the B-737 
aircraft in landing approach, t ouchdo^vn, rollout and turnoff for normal and CAT HE 
weather conditions. Preliminary results indicate that MLS can be used in place of 
ILS landing aids and that a single magnetic cable can be used for automated rollout 
and turnoff. Recommendations are made for further refinement of the model and 
additional testing to finalize a set of guidance laws for rollout and turnoff. 
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AmCRAFT DYNAMICS 


I, 

This section contains the derivation of the equations used in this study to 
model the aircraft motion during the landing approach covering the phases of 
capture, decrab and flare as well as the landing gear and tire dynamics for touch- 
down, rollout and turnoff. 

(a) Choice of the Coordinate Systems . It is customary to utilize linearized equa- 
tions of motion (Ref. [l] ) describing the deviation of the aircraft from a constant 
speed, fixed flight path angle, wings level equilibrium condition, in the body-stability 
axes system. While this is sufficient for the landing approach phase, it is inadequate 
for the rollout and turnoH due to the large deviation in air speed, flight path angle 
and the large changes in yaw during turnoff. In order to accommodate the broad 
spectrum of conditions it was decided to use the full non-linear equations of motion 
in a runway cartesian coordinate system to describe the linear acceleration of the 
aircraft and to use the body axis and the Euler angle of roll, pitch and yaw to des- 
cribe the angular acceleration equations. 

Since the full non-linear variation of the aerodynamic coefficients is small 
over the range of Mach numbers and aircraft altitudes in landing, and since the 
aerodynamic forces during rollout and turnoS become relatively unimportant when 
compared to the lauding gear and thrust forces and moments, it was decided to use 
a simplified form of constant lift, drag and moment coefficients, fitted to the B-737 
for the landing simulation. Numerical values of these aerodynamic coefficients 
were obtained from Langley Research Center and are given in Table 1. (Seepgs. 10-11) 

The orientation and sign convention for the runway inertial cartesian coordi- 
nate system is given in Figure 1. The aircraft position vector, R(x, y, z), is 
measured relative to the origin fixed on the center line of the runway with the posi- 
tive X axis pointing forward along the runway center line. The z axis is positive 
down, and the y axis is positive to the right in a right-handed orthogonal coordi- 
nate system. 
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Figure 1. Runway (Inertial) Coordinates and AircraE Position, Velocity Vectors (R, R) 



The laiit vectors a^, a^, anda^ are the tmit body axis vectors, lying along 
thex, y, and z body axes, respectively, expressed in the runway coordinate system 
in terms of the Euler angles 0, o* ). See Figs. 2a and 2b. The a^ axis is 
along the aircraft center line, positive, forward. ? 
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The a axis is the lateral aircraft axis positive, pointing along the right wing. 
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(b) The Aircraft Equation of Motion, The aircraft inertial velocity vector, , 

* » * • 

R(x, y, z) is given by the time rate of change of the position vector 


dt 


R =H 


( 2 ) 


The vector velocity of the aircraft relative to the atmosphere y(v , v , v ) 

■ JL A o 

is given by the aircraft inertial velocity vector and the winds, 

V = R - W + G -i-S / 


(3) 


where 






W is tlie constant wind vector, 

G is the gust vector , and 
S is the wind shear vector. 

The simulation of wind, shears and gust was supplied by the Langley Flight 
Instrument Laboratory contained in the Langley FLTFAST computer program. 

The aircraft linear acceleration is given by, 
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is the effective drag coefficient 
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C is the effective lift coefficient 
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is the effective side force coefficient 


C=C«jS+C 6 +C 6 
YR E Ysp sp 






The control surface deflections are defined In the notation. The thrust magni 

tude, k, and the gear forces (F , F , F ) will be derived later. 

X y z 

The aerodynamic angles are given by, 
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a^. V 


( 5 ) 


The Euler angular velocities of the body are given by the usual kinematic 
equations (Ref, [l]) 


6 = cos^q - sin0 r 

« 

0 = (r cos^ + qsin^)/cos 0 

, • 

^ = p + sin 6 0 


The angular acceleration equations of motion are given by, 
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where is the effective aerodynamic pitching moment coeffiOientg 
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C is the effective aerodynamic rolling moment coefficient ^ 
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and C is the effective aerodynamic yawing moment coefficient, 
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(Where p and r are the roll and yaw in the stability axis system, 
0 b 


p^ = cos O' p + sino: r 
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The term k a^ is the pitching moment due to thrust. 


In order to integrate Equations (7b) and (7c) it is necessary to eliminate the 
rm in Equation (7b) and the p term 
plished by inverting a 2x2 matrix as follows: 
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The twelve differential equations given b? Equations (2), (4), (16), (7a) and 
(8) are integrated to obtain the required motion of the aircraft. In the simulation 
a third order Range Kutta numerical integration scheme, using a time interval of 
i/20 second for landing approach and l/50 second for rollout and turnoff was 
found to produce accurate results. 

(c) Landing Dynamics. The landing gear forces and moments are generated by 
the ground reaction to the tire and strut deflection and the strut compression rates. 
In order to compute these deflections and compression rates, it is necessary to de- 
scribe the relative positions of the tire contact points of the main and nose wheels 
with reference to aircraft center of mass in the runway inertial system. 


The relative coordinates of the three extended, uncompressed gears of the 
B-737 are given in Refei’ence [2], and shown in Table 2, measured from the lead- 
ing edge of the mean aerodynamic chord at the aircraft center line in the plane of 
the MAC. The cartesian coordinates of the ith tire contact point are given by, 


x._ = X + G. G. c+x. 
iG i 




" z + z. +<py. - Bx. 
iG 11 1 


( 9 ) 


No gear load or moment can exist until 

( 10 ) 

Since the compi'ession struts are preloaded, no strut deformation can occur 
until the tire compression load is equal to the ith strut preload. The preload values 
are given in Reference [2] and illustrated in Figures 3, 4. The landing gear geometry 
is shown in Figure 5a. 

Thus, if 

O^ZiG^.0825 (11a) 
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TABLE 1 (cOEtinued) 
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TABLE 2 

E XTENDED GEAR LOCATIONS 


(Distances in feet from c.g.) 



X 

y 

z 

Gear No. 

Right Mam Gear 

-6.0 

8.58 

9.78 

1 

Left Main Gear 

-6.0 

-8.58 

9.78 

2 

Nose Gear 

28.3 

0 

8.77 
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Figure 3. Main Gear Strut Compression Force. 
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Figure 5(a). Landing Gear Geometry 
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the strut reaction is a liaear function of the tire compression and we have. 


F, = - F . 

1 pi o0825 


(lib) 


Once the preload has been exceeded the strut compression obeys the adiabatic com- 
pression law; that is, when 


,0825 

iG 1 


we have, 


F. = - F . 
I pi 


, z.-,0825 


1.4 


(lie) 


To obtain the strut damping force we first compute the strut compression 


rate z ^ . 
Gi 


^Gi^ (z ->-py^“qx^) 


(12) 


The damping force on the ith strut is given by, 




(13) 


The net strut force, acting normal to the ground, is given by. 


= F. + F„. 
Ni i Di 


(14) 


The tire drag force may now be computed from the tire friction coefficients, 

( 16 ) 

The modeling of the tire friction coefficient is tal<en from References [2 and 
3 ] . For dry runwa 5 ' conditions, 
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^ Polling 

For wet runway conditions, 

+ 100 
^ G 

Polling 10V-+200 

Lr 

Thus, for high landing speeds and wet runways, 

^^^BOLLING 

and for taxi speeds ^elow 20 kts,), and wet runways. 
The tire friction curves are shown in Fig, (5b). 


( 16 ) 
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(17a) 


(17b) 


Braldng for less than the maximum drag coefficient can be achieved by apply- 
ing partial braise pressure as a control variable. This is developed in the guidance 
and control sections of the report. Braking is not available on the nose wheel. 


The landing gear lateral forces and steering is taken from the theory deve- 
loped by W.B. Horne in Reference Cs]. We first compute the tire compression 
coefficient T}.: 


.5 F. 
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17 . > . 09 
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the tire compression is given by. 
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(19b) . 


= 7 vf. 7], 
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The tire cornering force coefficient is given, by, 
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The lateral tire cornering force is given 


F.^ . = K. (7. + or “ e.) 
Li 1 i i 


where e . is the ground track angle given by, 


f . = tan 
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y+r(x^^ cos (T~y.siucr) 


- X 


7. - is the steering angle, 

a is the aircraft yaw relative to the runway yaw; 

The landing gear forces, as given abov e are in tlie runway system, 
gear forces from the three gears are givf -i by . 
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(20a) \ 


(20 b) 
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The net 
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H. NAVIGATION AND FILTEKING 

This section contains the r*escription of the navigation system. The B-737 
is equipped with body mounted accelerometers and rate gyros. Figure 6 contains 
a block diagram of the functional relationships of the input IMU data, the MLS and 
radar altimeter measurements and the navigation and filtering functions that are 
required to feed information to the Autolaiid guidance and pilot displays. 

An optimum utilization of this equipment would require a Kalman type filter 
that would be capable of estimating MLS biases, winds, etc. 

Such a system is outlined in Reference [4 ] , and further study of such a 

J:. 

system is highly recommended. However, for the purposes of the study it was 
found sufficient to utilize a simpler fixed gain complementary filter technique for 
landing approach, rollout and turnoff. The navigation and complementary filter 
equations are given in this section. 

(a) Navigation . The body mounted accelerometers measure the specific foi-ce 
acting on the aircraft in body system coordinates. To transform the specific forc^ into 
the runway system it is i.. scessary to obtain the aircraft altitude angles with respect | 
to the runway. Tliese angles are supplied by the aircraft inertial platform unit. 

Let 

0 - roll 

0 = pitch (25) 

^ = yaw (measured from magnetic north) 

^ = lunway yaw ( measured from magnetic north) 
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Then, 


T^(0) = 
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COS0 

“Sin0 


sin0 

COS0 


(26) 
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Tg(CT) = 
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I "Sincr 




sincr 
cos cr 
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0- = 
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( 26 ) 

coat. 


The dead reckoning estimate of the aircraft acceleration in the runway sys- 
tem is given by. 


r 





The complementary filter utilises this acceleration to produce the estimated 
aircraft position and velocity by adding a correction proportional to the residual 
between the MLS measured aix'craft position in the runway system and the air- 
craft estimate of its position. j . 

(b) Complementary Filter, It is not surprising that with a good estimate of the 
aircraft attitude and acceleration it is possible to integrate the accelerations and ob- 
tain a good estimate of the aircraft position and velocity. Under such conditions, the 
main purpose of navigation aids such as ILS, MLS, etc., is to correct for errors 
in the UVIU due to gyro drift, accelerometer biases, etc. Such a system is outlined 
in Reference [43 and utilizes a square root versic of the Kalman filter with exponen- 
tial correlated process noise.. Such a filter possesses time varying gains and re- 
quires considerable skill in design, formulation and execution. The Kalman filter 
estimates corrections to the integrated IMU aircraft state by optimally minimizing 
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the covariance matrix of the errors in the aircraft state and the instrument biases. 

A complementary filter can be designed to provide a similar set of non- varying 
gains, without any relationship to the uncertainty in the aircraft state or the biases 
in the IMU and MLS measurements. By adjusting the complementary filter gains to 
pass the relative low frequency of the true aircraft acceleration It is possible to filter 
out the high frequency noise in the MLS measurements. The main criterion for such 
a filter is to provide a stable estimator, one that decays disturbances exponentially 
and lets the low frequency motions persist, A block diagram of the complementary 
filter is shown in Figure 7, 

Let one of the coordinates of the E{t) vector be x.(t). The differential 
equations of the complementary filter for the x. coordinate are given by, 


L 




( 28 ) 


dt iS' 3 ^ MLS. 

i 


where 


x.(t) = i component of R(t) 


It is plain that if the residual zero, the complementary 

filter will provide the dead reckoning solution of the aircraft motion in the runway, 
system. Moreover, the system is linear with constant coefficients (K , K and 
K > the filter gains), with non. homogeneous forcing fimctions due to the noisy MLS 

O 

measurements and the aircraft estimate of acceleration obtained from the gyros 
and accelerometers. By choosing the values of K . K and K , the solution can 

X . A o 

be designed to be stable and to filter out the high frequency noise. 

Since the MLS measiirements are available only In a discrete form at 
fixed intervals (see the next section) and the output of the accelerometers are 
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COMPLEMENTARY FMER FOR NAVIGATJQN.WITH 
BODY MOUNTED ACCELEROMETERS 



Figure 7, Complementary Filter for Navigation With Body Mounted Accelerometers Using MLS Measurements 









available at fixed intervals, we may produce a pseudo MLS measurement which is 
continuous and is equal to the actual MLS measurement at the discrete sampling 
time. 


We choose (t) to be: 

i 

(29) 


(30) 

i i 

With this definition of the MLS measurement we can effect an explicit solution of 
system of differential equations (Equations (28)) which has the required conditions 
that if the residual is zero we get the dead reckoning solution and which corrects the 
state to match the MLS observations, without the high frequency component. The 
solution is given by, 
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+ ( 0 i(jS^ - (j?) - J3(^^ + (})) sin (w AT)/w, 




2 2 „2 o 

a +ih +j8 -2o;;8 
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2a^j3 ^ cos (w AT)) 

+p? ) - {w^ + sin(^o AT)A) 




( 32 ) 


"3 = 


’K (O!(e““^'^»e“^^'^cos (wAT))' 

O 


— SAT 7 2 2 2 

.+ (oi-jS)e (sin (OJ AT)/co) / a +j3 +w -2^)3 


The constants, O!, jS, ^ are the roots of the characteristic equation of the 
system defined by Equations (28). 

(S + Oi)(S + iS+iW)(S-J- jS-iw)= S^+K_ S^+K.S+K = 0 

X A o 


“ 0J + 2jS 

Kg = 20!^ + + 

Kg = Oi (iS^ + 03^) 


(33) 


By choosing small positive values of a, jS, and to we can obtain a stable 
low pass filter. We recommend values of 


a= .08 


(34) 


OL 


03 ■-= 

^2 


To provide estimates of position, velocity and acceleration for the Autoland 
equations and display, we have. 


X,(t)=Xi2(t) 

X.(t) = X(t). + X.g(t) 


( 35 ) 
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not 


This section describes the equations in the B-737 Autoland guidance system 
that were used in the simulation study for landing approachc The section also derives 
a set of guidance equations for rollout and turnoff based on following a buried magne- 
tic cable. The rollout and turnoff guidance will operate using position and velocity 
information based on MLS measurements. Some question exists whether the range 
bias can be solved for and eliminated to provide for a safe high speed turn (30-50 
ft/sec) in CAT III conditions. 

(a) Landing Approach (Airborne) This section describes the guidance equations 
used in commanding a single banlced turn approach from an initial fixed azimuth and 
aircraft attitude to the runway azimuth at the point at which the aircraft reaches the 
runway centerline. The guidance consists of a vertical control law, a lateral contijol 
law, coupled with automatic throUle control to follow a circular ground track at a 
desired constant descent flight path angle. 

1) Vertical Path Control. The desired vertical descent rate is given 
by, 

HDOT = - sin (GAMMA) (36) 

C G 

where 

= 120 igsioTS / (36a) 

GAMMA = COMMAND PLIGHT PATH ANGLE 
The commanded vertical acceleration is zero, 

HDD^ = 0 (36b) 

The commanded altitude is 

(I) = H^ (I - 1) + H DOT^ * AT (36c) 

The vertical guidance equations are given by, 
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ROLL = 0 (DEGREES) 


QDEG = Q PEGREES/SEC) 

HDDOT = " Xg 
HDOT = 

EIGHT = X 

o 

== .16 pDOT^ - Y^), (Y (0) = HDOT ) 

[YjI ^2. 

A = .6 Y 

B = .09 (H (I) - EIGHT) 
c 

C == A+B+HDD 

c 

IcI*! 50. 

IF p (I) - EIGHT) <0 . f (a) - 1 . 
c 

D = Cf(Q!), (f (O') IS FIG. 8) 

E = -D+6. HDOT 

|e| ^5. 

Yg = .25 (E + HDDOT), (Y^ (0) ^ 0) 

IF STABILIZER TRIM IS OPERATIVE 
AND (E + HDDOT)^0., THEN 

Yg^O. 

Yg - . 0625 (QDEG -- Y^), (y^ (O) = 0) 

F = . 8 ((E + HDDOT) + Y + 2. 16 (QDEG - Y_) - . 004 ROLL^; 

o 3 

Y^ = 16. (F-Y^), (Y^(0) = 0) 


(37a) 


(37b) 

(37c) 

(37d) 
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( 38 ) 


The elevator command to the elevator servo is given by. 


DELEC = Y 



0 2 4 6 8 10 12 


a , Angle of Attack (deg) 

Figure 8, Vertical Path-Angle of Attack Function Control 

2) Horizontal Path Control, Let the desired ground path be a circular 

turn of radius of The aircraft is assumed to be approaching at a fixed azimuth 

of , The guidance law produces an aileron deflection designed to keep the air- 
o 

craft in a constant bank at a fixed airspeed of 120 knots. 

Let the aircx*aft ground track position be x , x and the turn circle center 

X ^ 

be Wj W^j measured in runway coordinates. The desired yav/ is given by, 

X y 

X “W 

* = tan" ~ — ^ (38a) 

5,-W 

2 y 

The cross track error is given by, 

CHTE = SIGN(Wy)(R^-y(5j^-W^)®-f(52-Wy)^ ) (36b) 
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The track angle error is given by, 


BEPEODUCIBU WY OP THE 
OBIGINAL I'Af;': i.i l.'A'A’i 


TANGE =- 




(3So) 


The desired bank angle is given in terms of the desired airspeed 



D 


tan 


-1 


“ 2 
V 
G 


G(SIGNCWy)R^) 


(38d) 


The guidance law bank command is given by. 


= & _C 
c D YD 


(X ^ +x/ )^TANGE - C CRTE 

i ^ JL 


where 


1^ I ^ 25.'^ 

' c ' 


^YD 




YD 


7,12 


(39) 


(39a) 


QPOT = .22 +120. 5/AIRSPEED 

The Autoland guidance equations for the banlc angle control is given below: 

Let 

A 

PDEG = P(DEGREES/SEC) 

ROLL = 0 (DEGREES) 

Y, = 40. (PDEG-Yg), (Yg (0) = 0) (40a) 

i^=5.(*/-Yg), (Yg(0) = 0) (40b) 

IF |5.(*^-Yg)| >4. 

f 

replace Equation (40b) by / 
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(40c) 


Yg=4. (SIGNCS^-Yg)) 

The aileron command to the aileron servo is given by. 


DEL AC= 1.82 QPOT (Y^ - - ROLL) 

a 0 

(41) 

|DELAC| ^ 

(41a) 


A test is made of when to set the desired bank angle to zero. This is done 
on the displacement from the runway center line, 

IF |x„ I ^ 100, FT 

BAmCD = 0, 

CRTS = X (42) 

A 

The guidance equations now return to the AUTO LAND capture mode for a 
conventional landing. This report does not contain the details of the capture, deerab, 
and flair logic. These modes were obtained from the Flight Division of Langley 
Research Center in the form of a computer program called FLTFAST.. The slm.ula'- 
tion of the Auto throttle thrust control, servo control of the stabilizer, elevator, 
aileron and speed brakes were talcen verbatim from the FLTFAST program de- 
signed to simulate the AUTOLAND control equations for the B-737. 

(b) Landing (Ground Guidance). The landing is divided into three phases, touch- 
down, rollout and turnoff. During the first phase, as soon as tiie main gear is pre- 
loaded, the throttle is automatically set in idle detent and the speed bralces are auto- 
matically engaged. Prior to touchdown this function is inhibited. 
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For the simulation, a two second linear slowout is employed during which 
time the ailerons, aileron spoilers, elevator, stabilizer and rudder are brought to 
zero and the control logic of decrab, flare, etc. is automatically- disconnected upon 
touchdown* The thrust control equations are also disengaged. 

For the two second period following touchdown, we have, , 

T^^^TIME^sTtj^ + 2. (43) 

DEmn(T )= RUDDER SERVO COMMAND 
^ AT TOUCHDOWN 

TIME - T 

XldC = 1. — (43a) 

The rudder servo command is given by a linear mix of DEIRC (TD) and the signal 
provided by the rollout guidance. 

DELR C = X KK (D EIR C (TD)) + (1 - XKK) (DELR C (ROLLOUT)) (43b) 

The guidance logic for the rudder control is given by a feedback of both the error in 
yaw and the error in yaw rate. 

Let the desired yaw, yaw rate, y deviation and lateral velocity be given by 
EDES = YAWD = YDES = VDES = 0. (44) 

Let the estimated values of these variables be 

A 

YAW = ^ (DEGREES) 

RDEG= R (DEGRESS/sec) . 

it f44a) 

YDES= Y, VDES = Y 

The rollout guidance follows; 


Y ^ = 10. (^EG "BDES"Y^) (Y^ (0) 0 . ) (44b) 

Yg= 2.5(Y-VDES-yg) (Y^ (0) = 0 ) (44c) 

ERUD = 7. Y + 4. (YAW -YAWD) + 2. (. l(Y-YDES) + Y ) (44d) 

IeRUh! ^10. ^ (44e) 
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EDAMP = RUDDER DAMPER COMMAND (FROM AUTOLAND EQU.) 
IedampI ^4. 

The rudder command to the servo during rollout is given by, 

DELRO (ROLLOUT) = ERUD + EDAMP (44f) 

The speed brakes flaps are brought to a full 40° deflection using a first 
order servo lag simulation 

Yg = 2. (40 - Y^) (45) 


The simulation of the thrust to the idle position (« 2000 „#) is given by a first 
order servo lag 

I^et THRUST (ID) be the desired thrust idle, then, 

-1 At 

TITRUST = e (THRUST-THRUST (ID)) + THRUST (ID) (46) 

(c) Rollout. Rollout begins at the end of the two second slowout. The thrust 
control is set to maximum reverse thrust (-17,400). Brakes are automatically 
applied by a law designed to bring the aircraft to a fixed taxi speed. Prior to reach- 
ing the taxi speed, when the aircraft reaches a ground speed of 30 knots, the thrust 
control is commanded to attain a thrust required to balance the drag force at the taxi 
speed and the rolling friction of the landing gear. 


Let 


VTAXI = desired taxi-speed 

An A n A 

THRUSTC(TAXI) = (VTAXI) + . 015 WEIGHT 


(47) 


Thus, when the aircraft speed reaches 30 knots, the thrust is given by, 

- 6AT 

THRUST -e * (THRUST-THRUST C (TAXI)) + THRUST C (TAXI) (48) 


The braking logic is given as follows: 
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lt52>KODUCn'n.lTY OF TTf 

ORIGINAL 1*^ FOv)i-l 


DLV = VTAXI-VG ^ 

. . . ^ (49a) 

Di VG = ‘’‘XgXg)/VG 

A CODES = (VG + VTAXI) DLV/150 - (DTVG - .015 G) 

Since the automatic sMd device controls the brake pressure to prevent skidding 
we may simulate the bralce pressure by modifying the effective coefficient of friction of 
the main gear. 


IF ACCDES ^ 0 ^Oa) 

DMUD = 0. (BRAICE COMMAND ZERO) 


IF ACCDES >0 

DMUD = ACCDES/G (50b) 

-2 At 

YMUD = e * (YMUD-DMUD)+DMUD 

|ymud| ^,5 

The coefficient of friction, simulating brake pressure, is now given by 


and 


Polling runways) 

^ VG + 100 YMUD 

^ ^ROLLING lOVG + 200 YMUD ^ runways) 


(50d) 

(50e) 


(d) Turnoff. The mechanism for providing steering information in this studj’’ is 
a single underground magnetic cable. The design and characteristics are taken 
from Reference fe]. In this section we derive the guidance law. Steering forces 
and moments are provided by the nose wheel as indicated in the earlier sections. 

For automatic control, the nose wheel is locked into a fixed turning ratio to the 
rudder servo as soon as the nose wheel vertical load exceeds the nose wheel pre~ 
load. At the high landing speeds when the rudder is still effective, the nose wheel 
forces are not significant. As the aircraft reaches taxi and turning speeds, the 
rudder loses its aerodynamic effectiveness and the nose wheel becomes the dominant 
steering mechanism. The ratio of the nose \vheel rotation to rudder rotation on the 
B-737 is given by, 
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( 51 ) 


ypEGREES) = -7~ 6^ pEGREES) 

Zb R 

The guidance law for turnoff uses the same g’ain constants as for the rollout. 
The estimates of the errors iu yaw and yaw rate are provided by different equations. 

A turn is a circular arc of a fixed radius. For the study the radius was 
taken to be 

R^ = 120. ft (52) 


As the aircraft approaches the start of the turnoff, a radio signal indicates the 
approach of the turn initiation and actuates the turn control circuit. 

Let the aircraft position at the start of the turn be x (fc ), then the center of 

Jm 5 

the turn circle is 


W 

X 


(t ) 

1 s' 


W =R^*SIGN 
y J- 


(53) 


The sign of the turn (to the right or left) is transmitted by a polarity code with the 
initiation signal. 

The magnetic cable yields a signal proportional to the error in yaw and the 
error in lateral deviation. In order to supply the information for the desired yaw 
rate and desired lateral deflection rate error, the information is generated from 
the assumption of a circular turn. 


For a turn of fixed radius, the desired yaw rate is given by 


RDES * (SIGN) (53a) 

To obtain a signal proportional to the track angle error we proceed as 
follows: 
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(53b) 


YS = Y “ SIGN * 

T 

DRW = ^ 

The error in deviation rate is given by 

• A -- 

(Y -VDES) = - (XSiXj^ + YSiXg)* (64) 

The turnoff guidance law follows the guidance logic given in Eq. (44a) to (44f), 


Reference [6] describes a steering law designed for the CV-880, which uses 
a position error system. Insufficient studies were carried out to determine the ad- 
vantage or disadvantage of the two laws. More studies are necessary before a firm 
recommendation can be made. 

Once the turn has been completed and the aircraft has a new heading, the 
guidance law returns to the straight line mode. The logic to initiate the new phase 
is a radio signal at the start of the straight line segment. The aircraft is given the 
X and X position and the yaw heading of the start of the center line of the taxi 

Jm ^ 

strip as a fixed code input. 


The desired yaw is now 


D RUNWAY 


(55a) 
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The desired yaw rate is 


RDSS = 0 , 


(55b) 


39 







hot YSlMM 


IV. THE MLS OBSERVABLES 

The characteristics of the MLS observables are described in documents 
supplied by Langley Research Center in Reference [7], The system used is for 
conical coordinates of asimuth (/S), elevation {a) and range. The range and aximuth 
origin is at the center of azimuth antenna facing outward toward the approaching air- 
craft along the runway centerline. The elevation origin is offset from the runway and 
is at the center of the elevation antenna. Thus, the cartesian coordinate system for 
the observable is the negative 01 the rum 'ay inertial system in the x and z 
direction and In the same direction as the y of the runway. We have, in the runway 
coordinate system, 


' X 


X 


'x' 

y 


-y 

+ 

-y 

z 


_z 


- ^ - 


MLS RUNWAY ANTENNA 


(56) 


. he origin of the MLS system is at the azimuth antenna and the position of the eleva- 
tion antenna is measured relative to it. See Figure 9(a). 


Let the MLS coordinates of the aircraft be X,,, Then the obser- 

M M * M 

vables in the MLS coordinate system are 




RANGE = J + Y,^^ + 

M M M 


)S= sin"^ (-Y^ /RANGE) 


0 ;= sin 


Z 


M 




(57) 


In the above, the z coordinate of the elevation antenna is assumed to be 

at the same altitude as the azimuth antenna, and and are the coordi- 

OE OE 

nates of the elevation antenna from the azimuth origin, in the MLS system. 
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AIRCRAFT 



GROUND 

PROJECTION 


Figure 9(a). The MLS Coordinate System. 
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To invert the solution and obtain X^, Y^, from range, ^ and a, 
we have 

Ym = -range (sin iS) (58) 

Zm = J (RANGE)" - 

where 

2 

g = Xog (sin O') 

Finally, to obtain the aircraft coordinates in the runway inertial system, we 

have, 


■ X 


X 


fx 1 

y 


y 

- 


1 

SUNWAY 

- 2 . 

AZIMUTH ANT 

- 


The noise characteristics are obtained from Reference [?] . 

For the exponential correlated noise ■■i iari:;o, and a, we have 

’’eF‘>=^el7^-V" 
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where u(t) is a random number of Kero mean, cr. is the standard deviation of 
the ith measurement and . 

_ At 

A. = e '^i (61a) 

For initialization 

n-(0) = u(0) (61b) 

Table 3 contains the values of the noise characteristics used in the MLS error 

model. 

To simulate data dropout, we use a random number with uniform distri- 
bution from -1 to 1 and set the observation equal to zero whenever the random 
number exceeds . 98. Thus, 2 percent of each observation string is lost. Since 
there are three observation types, 6% of the pseudo state observations are lost. 

This percentage can be varied. 

To simulate bad data, a similar random number with a uniform distribution 
is used, and an average percentage of the observations is set equal to a large num- 
ber (1000 times the bias) even though the date valid flag is set true. Thus bad 
data enters the system marked good a small percentage of the time. 

Biases are included in the simulation. Values of the biases are initially 
chosen from a random number generator and retained for the simulation run. 

The bias values are also given in Table 3. 

MLS "jitter" or "stagger" is simulated to produce the data train sequence 
as shown on Figure 9d, 

To accommodate data dropout, MLS invalids and bad data generally, the 
recommended course of action is to replace the bad data type by the aircraft esti- 
mate of the data missing (or considered rejectable) by the best estimate of the data 
type based on the current values of the vehicle state. Thus, if data is missing, we 

A A 

replace the bad data by range, azimuth or elevation from Eq. (57), using 
and 7 ,^^ computed from Eq. (56), This mode is still ut’uor study. 


RETRO: 
ORIGINAL • 
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TABLE 3 

MLS ERROR MODELS 


■ 1, — ^ . .1 , — , IM - 1 ■ 

Function 

Model 

A 

1 

0t{—) 

sec 

r 

u(n) 

Elevation 

y(n) = r y 1 - A^ u(n) + Ay(n-l) 

"0!T 

e 

19,100 

7, 01x10^^ Deg, 

IGRS 

Azimuth 

y(n) = tJi - A^ u(n) + Ay(n"l) 

-olT 

e 

0.971 

-3 

5,10x10 Deg. 

IGRS 

Range 

y(n) = r y 1 » A^ u (n) + Ay(n-1) 

CD 

I 

1.013 

21,1 Ft. 

IGRS 

, 


Note: IGRS = Independent Gaussian Random Sample 

Initialization is achieved by setting y(0) = u(0). 

See Figures 9(b) and 9(c) for illustrations of MLS jitter in azimuth and 
elevation. 

See Figures 10(a), 10 (b), 10(c) for illustrations of the cori’elated noise 
in MLS jitter in azimuth and elevation. 


















Elevation Jitter (Sec of Arc) 
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Figure 9(d). Time Sequence for Conical MLS. 
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Azimuth Noise (rad) 





V, 


IVIAGNETIC LEADER CABLE 


A single underground magnetic leader cable, as described in Reference [5] , 
is a centerline guide for prescribing the aircraft turn. The system is capable of 
generating a signal proportional to the distance from the cable centerline and a 
signal proportional to the yaw error from the cable centerline yaw. These are 
generated by functions of three voltages (Figure 11). 

V. 


CRTE=K, 


12 2 2 
V +V + V 
12 3 


(62a) 


TANGE = K, 


2 V. 


The sign of the heading can be obtained frona the phase between v and v . 

Z o 

The guidance law for the turn is similar to the rollout with the following 
changes ; 


RDES-RDEG= K, 


V V “V V 

3 2 2 3 


(62b) 


V 


yAWD-YAW = K, 


2 V, 


(62c) 


Finally, the rudder damper signal input to the rudder servo is shut off when the 
aircraft reaches the turn corner. The gain constants of the rudder -nose wheel 
guidance law are retained. 
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LEADER CABLE 



COIL 1 


Figure 11. Magnetic Leader Voltage Signals. 



VI. 


PRELIMINARY SIMULATION RESULTS 


This section contains a discussion of some of the preliminary results of the 
simulation program generated using the models outlined in this report. 

The airhorne guidance using MLS and the complementary filter successfully 
lands in the aircraft after executing a banked turn from a 165° yaw deviation from 
the runway with and without winds and gusts, Figures 12(a), 12(b), 12(c), 12(d) and 
12(e) are plots of a typical landing showing altitude, descent rate, pitch angle, run- 
way centerline offset and ground speed, versus time for capture, decrab, flair, 
touchdown, rollout and turnoff. 

The results indicate that the MLS system can be used for landing approach in 
mixed weather conditions. Additional studies, using the simulation program, should 
be carried out. The rollout and turnoff study shows that the aircraft can be automati- 
cally and rapidly decelerated to turnoff speeds using reverse thrust and braking. Tiie 
magnetic cable concept shows promise as a guidance tool aid. However, here also, 
additional studies must be carried out before concrete recommendations can be made. 

It is recommended that the simulation program described in this report be 
refined to Include a more accurate aerodynamic, landing gear and Autoland guidance 
model. Further, it is recommended that a noise model be included in the magnetic 
leader cable simulation. Finally, it is recommended that the competing rollout and 
turnoff guidance laws be studied so that a CAT III automated rollout and turnoff policy 
can be finalized for the B-737. 
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